Elusive electron-phonon coupling in quantitative analyses of the spectral function 
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We examine multiple techniques for extracting information from angle-resolved photoemission 
spectroscopy (ARPES) data, and test them against simulated spectral functions for electron-phonon 
coupling. We find that, in the low-coupling regime, it is possible to extract self-energy and bare- 
band parameters through a self-consistent Kramers-Kronig bare-band fitting routine and verify the 
momentum independence of the self-energy along the quasiparticle dispersion. We also show that 
the effective coupling parameters deduced from the renormalization of quasiparticle mass, velocity, 
and spectral weight are momentum dependent and, in general, distinct from the true microscopic 
coupling; the latter is thus not readily accessible in the quasiparticle dispersion revealed by ARPES. 
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Angle-resolved photoemission spectroscopy (ARPES) 
is a tool which, experimental difficulties aside, provides 
access to the electron-removal part of the momentum- 
resolved spectral function A(k, oj) [T]. This quantity is an 
extremely rich source of information since it depends on 
both the single-particle electronic dispersion e£ (the so- 
called 'bare-band') as well as the quasiparticle self-energy 
S(k, uj) = S'(k, bS) +iS"(k, u>), whose real and imaginary 
parts account for the energy renormalization and lifetime 
of an electron in a many-body system. The single-particle 
spectral function is generally written in the form: 

1 E"(k,w) 
( ' W) ~ ~n [to - el - E'(k, lu)Y + [E"(k, w)] 2 ' (1) 

In the case of k-independent self-energies, methods based 
on the Lorentzian fit of momentum distribution curves 
[MDCs, i.e. constant energy cuts of A(k, uS)\ have been 
used to extract £'(w) and E"(w) from the spectral func- 
tion measured by ARPES, and eventually to infer infor- 
mation on the nature and strength of the interactions 
dressing the quasiparticles in a variety of complex sys- 
tems P~HB]. However, the question of the validity of these 
approaches is still a pressing one because most methods 
hinge on some assumption and/or approximation for the 
bare-band e^. More fundamental, for instance in the case 
of electron-boson coupling, the commonly assumed link 
between coupling strength and quasiparticle renormal- 
ization, through the so-called mass enhancement factor 
(1 + A), is at best merely phenomenological and its gen- 
eral validity needs verification. This calls for a method- 
ological study based on a well-behaved and momentum- 
independent model self-energy. 

Here we investigate the possibility of extracting 
momentum-independent self-energies from A(k.,ui), with- 
out any a priori knowledge of the bare-band e£. We 
will test the performance of our approach on the spec- 
tral function generated with the high-order momentum- 
average approximation MA*"' [3 H] , for the single Hol- 
stein polaron model [S]: momentum- independent cou- 



pling between an optical phonon and a filled one-band 
system of non-interacting electrons. This is a highly over- 
simplified approach, in which strong interactions are not 
included as opposed to e.g. Ref.lTUl and even the effect 
of the Fermi sea is not accounted for (as there is no well- 
defined chemical potential, the latter will be positioned 
at the top of the first electron-removal state). Although 
not applicable to correlated electron systems, the plain 
Holstein polaron [S] is chosen as the minimalistic model 
to study the electron-phonon coupling problem. Since 
the MA( n ) has been shown to be extremely accurate ev- 
erywhere in parameter space [§] , it will allows us to study 
AQt, w) over a broad range of electron-phonon coupling. 

Before delving into the detailed self-energy analysis, 
let us illustrate the model Hamiltonian and emphasize 
some general aspects, relevant to the phenomenologi- 
cal description of spectral functions in terms of effective 
coupling and renormalization parameters. Limiting our- 
selves to the one dimensional case for simplicity (higher 
dimensions were found not to change the results qualita- 
tively), we have used the MA^") approximation to obtain 
self-consistent A(k, u) with highly accurate momentum- 
independent £(a>), from the Holstein Hamiltonian: 

^=E^ c ^+«E 6 q 6 q+^E c U^(^+6_ q ). 

k Q ViV fc,Q 

(2) 

Its terms describe, in order, an electron with disper- 
sion e\ = —2tcos(ka), an optical phonon with energy f2 
and momentum Q, and the on-site electron-phonon linear 
coupling [for N sites with periodic boundary conditions; 
c\ (c/c) and bg (6q) are the usual electron and phonon 
creation (annihilation) operators]. This leads to a dimen- 
sionless effective coupling X = g 2 /2tn. For this paper we 
set a = H = 1 and t = 50 meV, such that the bandwidth is 
200 meV and the Brillouin zone is 27rA wide; also note 
that an additional constant 1 meV FWHM Lorentzian 
broadening is applied, similar to an impurity scattering, 
to allow resolving the sharpest features in A(k,u>). 
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FIG. 1: (Color online), (a) A(k,uj) calculated within MA' 1 ' for f2 = 50meV and A = 0.5; the quasiparticle dispersion e\ and 
the bare-band e\ are also shown, (b) Quasiparticle and bare-band velocities, v k and v\, and (c) corresponding inverse masses, 
1/m^ and l/m\, according to the definitions v k — det/dk and l/m^ = d 2 e k /dk 2 . (d) Momentum-dependent quasiparticle 
renormalization as obtained from v^/v k , m q ,/m\, and the inverse quasiparticle coherence 1/ZJjf, where = f q A(k,Lu)doj is 
the quasiparticle-only integrated spectral weight; in the inset, these quantities are compared near fc = to the renormalization 
factors fl/W 9 and (1 + A), obtained from quasiparticle bandwidth W 9 and dimensionless coupling \ = g 2 /2tfl in our model. 



The spectral function calculated with MA*- 1 ' for ft = 
50meV and A = 0.5 is presented as a false-color plot in 
Fig.JTJa); as one can see, it deviates remarkably from 
the one of the uncoupled bare-band e\. The electron- 
removal spectrum is now comprised of a polaron quasi- 
particle band e k (i.e., the lowest-energy bound-state of 
the electron-phonon coupled system), and a continuum 
of excitations starting at an energy f2 below the top of 
the polaron band and roughly following the original lo- 
cation of the bare-band e\. Based on these simulated 
data and on the precise knowledge of the input value 
of A, it is possible to gauge the appropriateness of esti- 
mating the electron-phonon coupling strength from the 
observed renormalization of band velocities, masses, and 
quasiparticle coherence, as often done in the interpreta- 
tion of ARPES results Q]. To this end, in Fig.[ljb,c) we 
present the quasiparticle and bare-band velocities, v k and 
v k , and the inverse masses, l/m q k and {/m k (see caption 
of Fig. 1 for definitions). The corresponding ratios v k /v k 
and w^/rn^. , as well as the inverse quasiparticle coherence 
l/Z k (see again caption of Fig. [T]), are progressively larger 
than one the stronger the coupling strength; these quan- 
tities are usually equated to the renormalization or mass- 
enhancement factor (1+A), providing a potential path to 
the quantitative estimation of the electron-phonon cou- 
pling strength QUOD- 

As shown in Fig.JIJd), the velocity, mass, and spec- 
tral weight quasiparticle renormalizations are functions 
of k. While in general one could expect these quantities 
to all be distinct, we find that v k /v k = 1/Zf at all k, 
which is a direct consequence of the fc-independence of 
the self-energy derived from the Holstein Hamiltonian: 
by Taylor expanding £(fc, uj) in the vicinity of the quasi- 
particle pole, i.e. uj = e k + Suj, one can approximate the 
Green's function as G(k, uj) n Z^j(uj— e k ) in terms of the 



quasiparticle coherence Z k = 1/[1 — dT,(k,uj)/duj l^?], 
and obtain for the quasiparticle velocity = Z k [v\ + 
dT,(k,uj)/dk ]; the latter reduces to v k = Z k v k if 
<9£/9fc = 0. For the mass renormalization, we obtain 
m q k /m b k — 1/Z k only for k — and ir, which is simply 
a consequence of the fact that at the band extrema the 
velocities are zero and the corresponding rate of change 
away from the extrema has to follow the acceleration [see 
Fig.JIJd) and its inset; the horizontal dashed line empha- 
sizes the divergence of m q k /m k due to the vanishing of 
l/m q k \. Most importantly, v h k jv\ 1 m\/m b kl and 1/Z% can- 
not be compared directly to the momentum-independent 
renormalization factors fl/W q and (1+A), obtained from 
the quasiparticle bandwidth W q and the dimensionless 
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FIG. 2: (Color online). Renormalization parameters defined 
as in Fig.)!] plotted vs. the dimensionless coupling X = g 2 /2tQ. 
Note that the noise in v and 1/Z at k = ki originates from the 
determination of the inflection point fc;. 



coupling \ = g 2 /2tfl in our model [Fig.JTJd), inset]. Al- 
though coincidences in values can be observed at some 
electron momenta, these are purely circumstantial and 
cannot be generalized. 

The results in Fig. [I] demonstrate that in general one 
cannot directly extract a quantitative value for the cou- 
pling strength A from the observed renormalization of 
quasiparticle velocity, mass, and coherence (and this 
without even considering further complications originat- 
ing from the electron-correlation driven renormalization 
of electron-phonon coupling in real materials 1 1 3] ) . If 
one was to do that, the extracted value would vary with 
the chosen observable and the specific k value, making 
any conclusion arbitrary (e.g., at k — ir these renormaliza- 
tions overestimate A by a factor of almost 100). One may 
wonder, however, if any of these quantities scaled as (1+A) 
upon increasing A, so that if not the exact value at least 
the trend of the coupling strength could be captured, for 
example in an experiment performed as a function of dop- 
ing. In Fig. [2] we follow each of these quantities and the 
quasiparticle bandwidth as a function of A, contrasted 
against (1 + A). For momentum dependent quantities we 
must choose a k value: we plot Vq/vq, TOq/toq, and 1/Zq 
at k = where they all coincide, as well as v\./vl. and 
l/Z% at the inflection point k = fej of the quasiparticle 
band e| where rnf, /rn\, diverges. All quantities deviate 
dramatically from (1 + A) at large coupling values (e.g., 
overestimating the microscopic coupling by a factor of 
4 to 8 at A = 2), and are also rather poor indicators of 
the coupling strength in the low-coupling regime (A < 1; 
Fig. [2] inset). Interestingly, the inflection point velocity 
renormalization v\.Jv\. is linear in A over a range that 
scales with the ratio 0,/t. However, even this term devi- 
ates from (1+A) in a way dependent on the details of the 
model, such as the shape of the bare-band e\ ; in general, 
it cannot be used for the quantitative estimate of A. 

Determining e\ is a key step also in the attempt of ex- 
tracting real and imaginary parts of the self-energy from 
A(k,u), which will be the focus of the remainder of the 
paper. We analyze A(k,ui) in terms of MDCs at constant 
energy Li = Q. Since the self-energy in the present model 
is fc-independent, i.e. E^ and are constant, as long as 
e b k can be linearized in the vicinity of the MDC peak max- 
imum at k = k m , the MDC lineshape is Lorentzian. Note 
that the converse is not true: a Lorentzian MDC line- 
shape is not a sufficient condition to conclude £ = E(w) 
[13] , although the overlap of v\jv\ and l/Z k (Fig.[l]) con- 
firms the fc-independence along the quasiparticle disper- 
sion. By Taylor expanding e\ about the MDC peak max- 
imum at k = k m , i.e. e\ =s\ m + v \ m ■ {k — k m ) + ..., and 
noticing that a> = £j, m +E~, we can rewrite Eq.[l]as: 

Mk) -^(k-k m y + (Ak m r (3) 
where Ak m = — E^/uJj is the half- width half-maximum 
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FIG. 3: (Color online), (a) A(k,u) calculated within MA (1) 
for fl — 50meV and A = 0.1; also shown are the km path of 
MDC maxima, as well as the known bare-band and the one 
found through the KKBF analysis. (b,c) Real and imaginary 
part of the self-energy from the model (Ejj. ntmm ), the bare- 
band and MDC fitting routine (Ea/dc), and the KK trans- 
form of Xl D c (E' KK ) and Z' MDC (E£ K ). In (c) the MDC 
ratio results, E" at4o = — Ak m /Ao, are also shown. 

of the Lorentzian MDC and A = l/v\ = jA Q (k)dk. 
If the bare-band is not known it is possible to fit it, 
within an arbitrary offset, to any functional form which 
provides a value and derivative using a Kramers-Kronig 
bare-band fitting (KKBF) routine. This is done by first 
tracking k m and Ak m for every Co through a Lorentzian 
fit of the MDC (Eq. [3]), and then choosing e b k parameters 
such that Y,' MDC = uj — e b km and E^ /£)C = — Ak m are 
self-consistent with T,' KK and ^" KK calculated from the 
Kramers-Kronig (KK) relations: 

E^(fc, W ) = ±-v r o ^'mdcm (4) 

7T J_ 00 £-U) 
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FIG. 4: (Color online), (a) Various estimates for the imagi- 
nary part of the self-energy, defined as in the caption of Fig.[3] 
for A(k,u>) calculated within MA 1 - 1 ' for = 50meV and 
A = 0.5. (b) Deviation (i.e., average of the squared difference 
at each k m ) between estimated and known self-energies vs. A. 



In our implementation (Fig. [3]), a simple third order 
polynomial was used with an initial guess found by fit- 
ting MDC peak maxima. We then used the Levenberg- 
Marquardt Algorithm as implemented in the mpfit pack- 
age for IDL to vary band parameters. We found that 
the standard sum-of-squares minimization did not per- 
form as well as a concave-down function, since it placed 
too much weight on outlying points far away. In order to 
evaluate the integrals in Eq.|4]with a finite region of data, 
biased inverse polynomial fits where used to extrapolate 
tails before a Fourier-based transform was performed. 

The method outlined here varies slightly from tech- 
niques previously described in the literature, which of- 
ten deal with data very close to the Fermi energy and 
generally reduce the possible functional forms for e b k sub- 
stantially [T5J US] . While these methods present an exact 
solution for A(k, ui) based on a reduced bare-band, the 
present approach imposes no restrictions on e b k (other 
than it be differentiable near k m ), allowing fitting based 
on a wider variety of bare-band models and using KKBF 
to vary e b k parameters. The results of the KKBF proce- 
dure are presented for A = 0.1 in Fig. [3] and show that, in 
the low-coupling regime, the found bare-band and self- 
energies (Emdc) agree well with the known quantities 
from our model (as long as the MDC analysis is applica- 
ble: k m is far from band extrema where velocities vanish, 
and MDC widths are suitably small that the bare-band 
may be approximated as linear around k m ). Note that 
the self-energies are evaluated along the (fc m ,w)-path of 
the MDC maxima, along which the analysis is performed; 
this path deviates significantly from either e k or e\ close 
to the sharper one- and two-phonon structures [Fig.[3]ja)]. 
The reliability of the results is confirmed by the agree- 



ment between T,mdc and T.kk over the whole range, 
which provides an internal self-consistency check. 

Although very satisfactory at low coupling, at larger 
coupling the results from the KKBF procedure become 
progressively less accurate. This can be seen already for 
A = 0.5 in Fig.|4^a), where Y1' MDC and ^" KK fail to re- 
produce ^ known in the proximity of the sharp phonon- 
induced structures (S'kk nas & l so picked up different off- 
sets in the different flatter parts of the spectrum, making 
setting its overall offset difficult). For a more compre- 
hensive description of the deviation of these quantities 
from ^ known upon increasing A, in Fig.Qb) we present 
the average of the squared difference at each k m be- 
tween the estimated and known self-energies versus A: 
for both YI' MDC and YJ' KK , this estimate shows a rapid 
and consistent increase with A, indicating a failure of the 
KKBF analysis and similar methods in the intermediate 
to strong electron-phonon coupling regime. 

Before concluding, we will point out an alternative, 
possibly more practical approach, which allows us to 
tackle the problem over a larger range of A, at least 
for £". In our method, determining real and imaginary 
parts of the self-energy hinges on finding a proper expres- 
sion for the bare-band through the KKBF routine. The 



imaginary part however, 



Ak m , only requires 



knowledge of v\ , and this can be obtained directly from 
A(k, uj) in two independent ways. The first one is through 
the momentum integral of A^(k) in Eq.[3j which returns 
directly Aq = \/v\. ; this allows a simple estimate of S" 
[and equally accurate to KKBF's for A = 0.1, Fig.^c)] 
from the MDC width/integral ratio: S" atio = — Ak m /A . 
The second one is through the equivalence v\lv\ = \jZ\ 
discussed in the context of Fig.[ljd), and the possibility 
of estimating v k and Z k directly from the data as the 
momentum derivative and energy integral of the quasi- 
particle band e\ (Fig.[l]): E'^ DC = ~v q k Ak m /Z q k . The 
so-obtained E" atio and YI' MDC are compared in Fig.BTa); 
although especially E" otio deviates strongly from Ef„ 0!lm 
upon increasing A at the sharp phonon-induced features 
[Fig.Qb)], the general behavior is that when E"atio ~ 
Yl' MDC they also match ^ known almost exactly. Thus, 
S" Qtio and H'l lDC can be used to obtain a precise anchor 
mesh for E" over a range of energies, without the critical 
step of finding the bare-band through KKBF. 

In summary we have shown that, at variance with a 
common phenomenological practice in the interpretation 
of ARPES data, even in the simplest case of electron- 
phonon coupling described by the Holstein model it is not 
possible to obtain the microscopic coupling strength from 
the observed renormalization of quasiparticle coherence, 
velocity, mass, or bandwidth (e.g., v k ./vl. would overes- 
timate A by a factor of 2, 3, and 8 for A = 1, 1.5, and 
2). In this sense, the coupling A still remains an elusive 
quantity. Through the KKBF analysis we can however 
gain access to bare-band and self-energies, which if prop- 
erly modelled might provide information on the nature 
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and strength of the underlying interactions, at least in 
the momentum independent, low-coupling regime. 
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